Is the Galactic Cosmic Ray Spectrum Constant in Time? 
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ABSTRACT 

The hypothesis is considered that the present Galactic cosmic ray spectrum is 
at present softer than its time average due to source intermittency. Measurements 
of muogenic nuclides underground could provide an independent measurement of 
the time averaged spectrum. Source intermittency could also account for the 
surprising low anisotropy reported by the IceCube collaboration. Predictions 
for Galactic emission of ultrahigh-energy quanta, such as UHE gamma rays and 
neutrinos, might be higher or lower than previously estimated. 



1. Introduction 

The measured residence spectrum of Galactic cosmic rays (CR) is N r (E)dE oc E~ 2 ' 7 dE 
from 1 GeV to E^ nee = 3 PeV, and E~ 3 - 3 from E^ nee to E an ki e = 4 EeV. This is considerably 
steeper than the theoretical expectation for the injection spectrum of Ni(E)dE oc E~ 2a dE for 
diffusive shock acceleration (DSA) by strong shocks. Much of the difference is attributed to 
energy dependent escape rate R(E) from the Galaxy R(E) oc E d . However, it is questionable 
whether the entire difference is due to this energy dependence. For one thing, the secondary 
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spectrum is about N s (E)dE oc E~ 3 -°dE suggesting that d ~ 0.3 (Trotta et al., 2011), and, 
by implication, Ni(E)dE oc E~ 2A dE. Moreover, if d were much more than 0.3 over many 
decades of energy, then the anisotropy at the highest energy of Galactic CR near 1 EeV, 
would be much higher than the reported upper bound of of ~ 10~ 2 . 

There could be many reasons for the injection index to be much steeper than 2.0, but 
it would in any case be good to know its exact value, as the total energy requirements 
for Galactic CR are a sensitive function of it. Moreover, the UHE emission that could be 
expected from dense environments that offer a thick target to such CR are an extremely 
sensitive function of the injection index. In this letter we consider the unconventional pos- 
sibility that the time average residence spectrum is actually flatter than -2.7 by some small 
amount e. We would then conclude that the injection index is -(2.4- 2e). 

The possibility and consequences of time variability of CR have received increased at- 
tention (Melott and Thomas, 2011, Erlykhin and Wolfendale, 2010 and references within) 
following the suggestions that CR flux correlates with cloud cover (Svensmark, Bondo and 
Svensmark, 1997) and that very high energy CR trigger lightning (Gurevich, Milikh, and 
Roussel-Dupre 1992). Erlykhin and Wolfendale have specifically calculated the variation in 
CR intensity due to discrete, local supernovae at various energies. As far as we know, there 
is circumstantial evidence that cloud cover, lightning, and cosmic ray flux all anti-correlate 
with solar activity. There is not yet a firm causal link established between cosmic ray flux 
and lightning or cloud cover. The possibility remains that they could all be tied to a fourth 
correlate, however, it would still be a mystery what this fourth correlate would be. In any 
case, variations due to terrestrial and heliospheric conditions are separate from variations in 
the local Galactic flux. This paper is concerned with the latter. 

At energies much higher than 1 GeV, CR production in the Galaxy may be intermittent. 
Intermittency can lower the anisotropy relative to steady state during the lulls between the 
intermittent production episodes (Pohl and Eichler, 2011, Pohl and Eichler, 2012, Kumar 
and Eichler, 2012). The recently reported anisotropics and/or upper bounds may be at least 
partly understandable in these terms. This explanation for the unexpectedly low anisotropy 
would then predict a slightly steeper present day spectrum than the time averaged spectrum. 

Abundances of terrestrial cosmogenic nuclides (TCN) can set limits on recent time 
averaged variations of the cosmic ray flux at Earth, though not on short term variability. 
They demonstrate that the cosmic ray number flux in the Solar System could not have varied 
on average by more than about 20 or 30 percent of its presently measured value. (Some weak 
constraint on even short term cosmic ray variability is placed by the continued existence of 
life on Earth for hundreds of millions of years. The possibility remains that astrophysically 
brief fluctuations had interesting terrestrial impact that fell short of total lethality.) 
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The discreteness of supernovae, assuming they are the main source of CR, should in any 
case provide an a priori case for some time variability in the CR flux at Earth. Supernovae 
are distributed in the plane of the Galaxy, and most local cosmic rays come from sources 
at a distance commensurate with the halo size, H, conventionally taken to be a few kpc. 
Beyond H, cosmic rays preferentially diffuse out of the face of the Galaxy rather than across 
the disk, and their contribution at large distances, r, is exponentially small in r (Pohl and 
Eichler, 2011). 

Fluctuations around the time average should increase as the escape time from the Galaxy 
decreases. If the escape time decreases with energy, E, so then does N(E), and the level of 
fluctuation increases with E. This amounts to fluctuations in the spectral index. Taking for 
example the rate of supernovae, P, the CR source strength, Q, and diffusion with diffusion 
coefficient D, the average CR flux in the Galaxy (of effective radius R) is 

n (e) ~ ^ PTcsc m 

where r esc = H 2 /2D is the escape time. A single source provides CRs that are, let us assume 
for simplicity, are distributed over a sphere of radius r(t) = y/~Dt, where t is the time since 
CR release. For r(t) < H, the density of these new CRs therefore is 

N(E) 3Q = — (2) 

sl j -4vrr(t)3 47r(Z>t) 3/2 ' { ) 

If the supernova is at distance r, one expects a spike in CR flux of relative amplitude 

3R 2 H 

V ~ N ~4r 3 Pr ' [ ) 

Jv av ^" 1 'esc 

For GeV-band CRs P r csc 10 6 , and therefore we have to be within 

r c {rj) ~ (100 pc) ?T 1/3 \ D[l G(N) ) ( 4 ) 

of the source to expect a spike of amplitude r\. Note that r c is very weakly dependent on 
CR energy. The duration of a large fluctuation (77 > 1) CR event is about 

T - - 2r csc 4 - (6 • 10 4 yr) ( ^^Y^ ■ (5) 



D(E) H 2 v J ' \D(1 GeV) J 

The frequency of such CR events is 

2 / n(n\ \ 2 / 3 



f~pH~io-*p( D ^ ) (6) 

1 ~ R 2 ~ \r]D(l GeV) / 
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and the fraction of time, 77, that one can expect to experience an event is given by the product 
of frequency and duration, 



Thus, significant events with r] > 1 are somewhat exceptional, but not unexpected. They 
are likely to be accompanied by spectral distortions because, at a fixed distance from the 
source, high-energy CRs arrive earlier than low-energy CRs. 

In this letter, we consider the possibility that the average residence index in the Galaxy 
is -(2.7 - e), e ~ lO^ 1 , that the escape rate exponent S is in fact 0.3 + e, and that the time 
average CR injection index is -(2.4-2e) rather than -2.4. This time averaged spectrum is 
slightly "propped up" by nearby supernovae, which introduce freshly accelerated populations. 
That the present day spectrum is then as steep as E~ 2 - 7 , according to the present hypothesis, 
due to the paucity of recent, nearby supernovae. 



We suggest that the variation of the CR spectral index over time in the GeV to TeV 
range could affect, at detectable levels, the rate of muogenic nuclei production in deep rock 
formations that exceed 3 Myr in age. Studying nuclear cosmogenesis at and below the 
Earth's surface has the advantage, in searching for evidence of hard transient sources, that 
the overlying material filters out much more of the contribution from low energy primaries 
than from the high energy ones. In particular, muons typically arrive from CR primaries of 
0.03 < E ~ 1 TeV, and would be more prevalent, relative to the total cosmic ray flux, for 
harder spectra. 

The depth would need to be large enough to be unaffected by possible uncertainties in 
the erosion rate at the surface. Variations in the CR intensity (due to nearby supernovae, 
or due to long term magnetic field change in the sun or Earth) within the past 3 Myr 
could possibly be tested by 10 Be abundance, to measure the exposure age - and, at the same 
geographical location, at depths of at least several meters, to isolate the muogenic component. 
Alternatively, measurements deep enough to give reliable exposure ages, independent of 
uncertainties in erosion rates, would give an absolute measurement of the past muon flux. 
The cross section for muogenic 53 Mn (half life 3.7 Myr, exponential decay time 5.3 Myr.) 
from an iron target is 3.8 mb at = 190 GeV. (Heisinger et al, 2001). The production rate 
per target atom, with a muon flux that is the same as at sea level, is then 10~ 21 /yr (equation 
17 of Heisinger, et. al, 2001). Here we have allowed for the fact that the mean energy at 





2. Cosmogenic Nuclei Diagnostics 
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sea level is only ~ 10 GeV, and that the cross section for muogenic 53 Mn production varies 
as ~ E®' 86 . The production flux remains within an order of magnitude of the sea-level value 
down to depths of 10 4 g/cm 2 , and analysis of the 53 Mn content of subterranean iron ore would 
give a straightforward measurement of the exponentially-weighted average muon flux over 
the past 5 Myr. Comparison with muogenic 10 Be (half life 1.6 Myr, production cross section 
on O at 190 GeV of 0.094 mb), could help date the contribution of recent nearby supernovae 
within the residence time in the disk (~ 3 Myr). Evidence for a recent nearby supernova, 
occurring about 2.8 Myr ago, has been presented by Knie et al.(2004). 

Measuring changes in the CR spectrum that are restricted to energy range E p 3> 1 TeV 
would be more difficult, so a time averaged CR flux that is well above the present only at 
E p >l TeV (such as might result from CR injection by ultrarelativistic shocks with Lorentz 
factors exceeding 100) would be hard to rule out. Muons from primaries at E p 3> 1 TeV 
dominate the total flux only at depths of > 10 5 g cm -2 , where the nuclear muogenesis rate 
per target atom is only about 10 _23 /yr. A net collection of 10 6 53 Mn atoms, accumulating 
over 10 6 5 yr, would thus require a purified sample of 10 22 5 iron nuclei, i.e. several grams. 

The dependence of muon secondary flux, and its depth dependence, can be estimated 
with a simple analytic model: Muons at sea level come mainly from primaries at E > 
lOOGeV, which create center-of-mass Lorentz factors that are larger than the ratio of the 
muon flight time to its rest frame lifetime. A detailed account of the physics is given in 
Dorman (2004), and a look-up table for muon production as a function of primary energy 
is given by Atri and Melott (2011). ; for the following discussion, a rough analytic estimate 
will suffice, and we take the atmosphere to be isothermal and planar: As primaries create 
the first generation muons at about 1 interaction depth which we take to be 85 g/ cm 2 in air, 
or about 1/12 of the total vertical grammage, 1030 g/cm 2 , of the atmosphere, this typically 
occurs at an altitude of h* ~ 18 /cos9 km, where 9 is the angle of incidence of the primary 
relative to the zenith. From the muon production altitude h* to altitude h, the proper flight 
time in the muon rest frame is given by 

ph ph 

r(9,h)= dr= -dti/~f{ti,9)c cos 9 (8) 



where 

rh 

7 (/i, 9) = 7 (/i„) - 20G(h, 9)/G(0, 0) = 7 (/i„) - 20 ' 



p(h')dh'/cos9 



h 



/G(0,0) (9) 



The above expression assumes that muons lose 20 times their rest energy, 1.05 • 10 8 eV, 
per vertical atmospheric grammage G(0,0) - i.e. about 2 MeV per gram over 1.03 • 10 3 g 
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cm 2 . Assuming a density distribution of p(h) = p exp[—h/h s ], neglecting terms of order 
exp[—h*/h a ], and integrating equation (jSJ) from h* to altitude h, h <C h*, we obtain 

fa h\ h ^ c i [e fe -^ 7 ,cosg-20] 
^' = 7^0^ [e*/*. 7 . cos 0-20] 
where 7* = 7(/i*). For /i s ~ 8 km, a reasonable value for the Earth's atmosphere over the 
lowest 3 scale heights, we can write h s /c ~ 12^, where r M is the muon lifetime. The decay 
factor is then 

-((e^-20/ 7 *cos^ 12/7 * cose 



e -T(h,e)/T„ 



( e h«/h s _ 20/7, cos e 

The muon flux at sea level, ^(7(0), 9, 0), can be expressed as 



11) 



Integrating over \i and over 7*, over the range where 7* > 12///, and assuming that the 



d 2 F^,^h«) ^ d 2 F fl (l,fi,h t ) ^ -2-p — 

sky muon flux at sea level, F M , 



muons have a spectrum of — M ^*^' — ~ — dd 2 * '* 2 p = C*7* 2 p , we obtain for the all 



p 1 /»oo 

F M (0)~C / dcostf / 7*^ 

JO 



on / /1\ "1 12/7» cos 

20/7 - COS(,) I (13) 



( e h/h s _ 20 / 7* cos $) 

where the lower limit of the integral is determined by the condition that the integrand be 
positive. 

When 7* ^> h*/ cr^ cos 9 ~ 27/ cos 9, r(h, 6*)/r M <C 1, and the decay factor is about unity. 
Making the approximation that the decay factor is close to when 7* < h*/h s ~ 27/ cos 9, 
and 1 when 7* > 27/ cos#, we can estimate the above integral as 

f 1 /*oo 

F M (0)~C/ dcosfl / 7" 2 "^7* = [l/(p + l)(p + 2)(27) p+1 ]C (14) 

JO J27/cos6» 

Thus, decreasing p from 0.7 to 0.6 with ■^■(lGreV) remaining constant, increases the 
flux near sea level, g = G{h, 0)/G(0, 0) ~ 1, by a factor of (27) ai (1.7/1.6)(2.7/2.6) ~ 1.53 . 
This is a good enough estimate for considering the general feasibility of using deep muogenic 
isotopes to measure the past cosmic ray spectrum. 

Deep underground, where g S> 27/20, the incident muons need to begin with a minimum 
7* of 20g/ cos 9 in order to survive energy loss, even though their decay is negligible. The 
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flux can then be approximated as 




20g/ cosd 



OO 



(15) 



A change in the spectral index -2-p from -2.7 to -2.6 then lead to a change in the all sky 
muon flux by a factor of (20#) ai (1.7/1.6)(2.7/2.6) ~ 1.5# 01 . 



We now consider how much the observed spectral index is likely to vary with time 
at Earth. We have calculated an expected intensity vs. time plot for randomly placed 
CR sources in the Galaxy, weighted according to location in proportion to local starlight 
production. The two figures plot the "effective" escape index (EEI) - i.e. the log of the 
ratio, as a function of energy, of the CR flux to the source spectrum at the same energy 
- assuming a value for d of 1/3 and assuming that the cosmic ray sources are dominated 
by bright supernovae. This index differential can vary from —0.2, within the spiral arm, to 
—0.45 between spiral arms. This is discussed in more detail elsewhere (Kumar and Eichler, 
in preparation). Over timescales of 3 Myr, the variation could be of the order of 0.05, which 
would produce a 53 Mn anomaly of > 10 percent. This variation would be detectable if the 
muon flux averaged over the 53 Mn lifetime could be measured to ~ 10 percent accuracy at 
g ~ 100. 

The solar system now lies between two major spiral arms but inside a minor spiral 
arm. There are certainly several nearby, recent supernovae, and they and others may have 
contributed to some fluctuation of the spectral index over the past several million years. 
In addition, however, the excursion into and out of major spiral arms, as illustrated in the 
figures, may introduce variations of as much as 0.1 or larger in the effective escape index 
around its average value. The timescale of this variation is longer than several million years, 
and would not be manifest in 53 Mn anomalies, but perhaps would be in longer lived isotopes. 
If the CR residence time in the halo is larger than the residence time of the Solar System 
inside a spiral arm, then the time-averaged value of the EEI is perhaps closer to -0.4 than to 
0.3, and the implied source spectral index is flatter by 0.2 than implied by an EEI of -0.3. 
That is, the inferred source spectrum in closer to -2.2 than -2.4. 



3. Predicted Variation in Spectral Index 
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time in Myr 



Fig. 1. — CRs flux variation at four different energies for isotropic diffusion in disk geometry, 
as sun crosses spiral arms at about 120 Myr (relative time) and 260 Myr. Diffusion coefficient 
for CRs of energy E is taken to be 6 x 10 27 (E/lGeV) 1 ^ 3 cm 2 s~ l and disk boundary is at 1 
Kpc. Top and bottom panels are for source rates 10 -2 and 10~ 3 yr _1 respectively. The blue, 
green, red, and cyan lines are for 1 GeV, 10 GeV, 100 GeV, and 1 TeV respectively, in each 
panel. 



- 9 - 




100 120 140 160 180 200 220 240 260 280 300 
time in Myr 



Fig. 2. — Temporal variation of difference in residence and injection indices, considering the 
CRs spectrum between 1 GeV and 1 TeV (figure [3]). As in the figure [31 top and bottom 
panels are for source rate 10 -2 and 10 -3 yr _1 respectively. Sub-panel in the top panel shows 
the variation between 120 Myr and 130 Myr (a period when Earth lies just inside a spiral 
arm) on a smaller time scale (0.1 Myr). 
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4. Further Discussion 

The expected emission by the Galaxy of ultrahigh-energy neutrinos, typically generated 
by primaries in excess of ~ 10 14 eV, would be raised by about a factor of 3 for each change 
of 0.1 in the EEI, and the contrast between the spiral arms and the spaces between them 
might be enhanced by an even larger factor. According to figure 1, this factor could be as 
high as an order of magnitude or more at such high energies. Similar considerations apply 
to UHE gamma ray emission from young sources. 

The hypothesis that we are living in an unusual era, when the CR flux at high energies is 
well below its average value in the Galaxy, might be motivated by anthropic considerations: 
e.g., if a normal CR flux were somehow detrimental to intelligent life. The present level of 
ionizing radiation from CR flux is less than the component due to terrestrial radon, so it is 
hard to see how changes of less than a factor of 2 in the CR flux could have serious astrobio- 
logical consequences, though it cannot be ruled out beyond reasonable doubt. At very high 
energy (~ 10 2 PeV), intermittency can cause much larger fluctuations in the CR intensity, 
and if high energy airshowers have astrobiological connections (such as lightning, rainfall and 
nitrate formation at low altitudes) then it is conceivable that the not unexpected changes 
in the spectral index of CR could conceivably affect terrestrial affairs and the development 
of life and civilization. Even if very high energy CR are produced by the same events as 
lower en energy CR, changes in the spectral index of order 0.1 are not unexpected, and even 
this modest change could influence the flux at 100 PeV by nearly an order of magnitude. A 
recent nearby supernova, for example, could have briefly raised the CR flux enough to have 
encouraged life (e.g. by encouraging lightning and, in turn amino acid and/or nitrate forma- 
tion), or destroyed it by climatic disturbance or excessive radiation and biological mutation 
rate. 

Neutrino emission from the spiral arms, and UHE gamma ray emission from supernovae 
remnants are sensitive functions of their respective local CR spectral index. For a given 
gamma ray flux at 10 2 MeV, they could set useful limits on the spatial variation of the CR 
spectral index. 
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